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To accelerate functional genomic research in tomato, we 
developed a Micro-Tom TILLING (Targeting Induced Local 
Lesions In Genomes) platform. DNA pools were con- 
structed from 3,052 ethyl methanesulfonate (EMS) 
mutant lines treated with 0.5 or 1.0% EMS. The mutation 
frequency was calculated by screening 10 genes. The 0.5% 
EMS population had a mild mutation frequency of one 
mutation per 1,710 kb, whereas the 1.0% EMS population 
had a frequency of one mutation per 737 kb, a frequency 
suitable for producing an allelic series of mutations in the 
target genes. The overall mutation frequency was one 
mutation per 1,237 kb, which affected an average of 
three alleles per kilobase screened. To assess whether a 
Micro-Tom TILLING platform could be used for efficient 
mutant isolation, six ethylene receptor genes in tomato 
(S/E7R7-S/E7R6) were screened. Two allelic mutants of 
S/E7RT (Sletrl'l and Sletrl-l) that resulted in reduced 
ethylene responses were identified, indicating that our 
Micro-Tom TILLING platform provides a powerful tool 
for the rapid detection of mutations in an EMS mutant 
library. This work provides a practical and publicly ac- 
cessible tool for the study of fruit biology and for obtain- 
ing novel genetic material that can be used to improve 
important agronomic traits in tomato. 

Keywords: Ethylene receptor • Fruit ripening • Micro-Tom • 
Reverse genetics • TILLING • Tomato. 

Abbreviations: CAPS, cleaved amplified polymorphic se- 
quence; CODDLE, codons optimized to discover deleterious 
lesions; DAF, days after flowering; EMS, ethyl methanesulfo- 
nate; EST, expressed sequence tag; GABA, y-aminobutyric 
acid; NBRP, National BioResource Project; SOL, The 
International Solanaceae Genomics Project; TILLING, 
Targeting Induced Local Lesions In Genomes. 



Introduction 

Tomato {Solanum lycopersicum) is a vegetable crop that is 
commercially valuable worldwide. It is a model system for 
studying many aspects of fruit biology, including develop- 
ment, ripening and metabolism (Giovannoni 2004, Tanksley 
2004, Carrari and Fernie 2006, Giovannoni 2007, Barry and 
Giovannoni 2007, Pineda et al. 2010). Insight into the carot- 
enoid biosynthetic pathway, chromoplast differentiation 
(Isaacson et al. 2002, Galpaz et al. 2008, Egea et al. 2010, 
Stigliani et al. 2011) and y-aminobutyric acid (GABA) metab- 
olism (Akihiro et al. 2008, Yin et al. 2010b) has been obtained 
using tomato as a model plant. Tomato has been widely used 
for studies of the biotic and abiotic response, such as salinity 
and wounding stress (Yin et al. 2010a, Sun et al. 2010, Sato 
et al. 2011) and infection by plant-parasitic nematodes 
(Bhattarai et al. 2008, Uehara et al. 2010). Large-scale expressed 
sequence tag (EST) analysis demonstrated that approximately 
30% of expressed genes in tomato had no significant similarity 
to Arabidopsis genes (Van der Hoeven et al. 2002). Tomato is a 
suitable system for studying unique biological phenomena not 
harbored by Arabidopsis. 

The tomato genome was entirely sequenced by The 
International Solanaceae Genomics Project (SOL), and many 
of the gene sequences can be retrieved from databases 
(Mueller et al. 2005a, Mueller et al. 2005b, Mueller et al. 
2009). Now, it is important to use the genome sequence infor- 
mation to elucidate the gene functions. TILLING (Targeting 
Induced Local Lesions In Genomes) allows for the identification 
of an allelic series of mutants with a range of modified func- 
tions for desired genes (McCallum et al. 2000a, McCallum et al. 
2000b, Colbert et al. 2001, Comai and Henikoff 2006). TILLING 
has been applied in many model and crop species, with no 
exception of tomato. In contrast to transgenic methods. 
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mutagenesis has the advantage of producing mutant alleles 
directly for use in breeding programs as genetic sources or 
markers. Ethyl methanesulfonate (EMS) and fast-neutron 
mutant populations have been created in the M82 genetic 
background. Mutants were classified according to their visible 
phenotype, and this information is publicly available through 
a database (Menda et al. 2004). Using this M82 EMS mutant 
population, a TILLING platform was recently created by an- 
other group (Piron et al. 2010). Other cultivars, such as cv. 
Red Setter and cv. TPAADASU (Gady et al. 2009, Minoia 
et al. 2010), also served as material for creating EMS popula- 
tions, and TILLING platforms were developed from these 
mutant populations. 

Recently, a dwarf tomato cv. Micro-Tom has been exten- 
sively used as a model variety to accelerate functional genomic 
research (lijima et al. 2008, Aoki et al. 2010). Like Arabidopsis, 
Micro-Tom is a suitable model plant due to its small plant 
size and short life cycle, and the existence of high efficiency 
transformation methods and indoor cultivation protocols 
under normal fluorescent conditions (Meissner et al. 1997, 
Emmanuel and Levy 2002, Sun et al. 2006, Watanabe et al. 
2007). To establish functional genomic tools in tomato, we 
previously developed nearly saturated Micro-Tom EMS and 
y-ray-irradiated mutant populations (Matsukura et al. 2007, 
Watanabe et al. 2007, Saito et al. 2011). Information on this 
mutant collection is available in the TOMATOMA database 
(http://tomatoma.nbrp.jp/). Although numerous tomato mu- 
tant populations have been generated by several groups, re- 
ports on large-scale mutant screening are limited, probably 
due to the large size of common tomato cultivars used in 
those studies. The cultivation of common cultivars requires 
extensive space in the greenhouse or field, making it difficult 
for most researchers to characterize even a single mutant. 
Therefore, a TILLING platform made from the Micro-Tom 
mutant population would greatly facilitate the functional ana- 
lysis of tomato genes. 

Several fruit-ripening mutants have been reported in 
tomato, such as ripening-inhibitor (rin), Colorless non-ripening 
(Cnr), non-ripening (nor), Green-ripe (Gr) and Nei/er-ripe (Nr). 
The causative genes were isolated and shown to act as up- 
stream regulators of the ethylene signaling network or in ethyl- 
ene perception (Wilkinson et al. 1995, Vrebalov et al. 2002, 
Giovannoni 2004, Barry et al. 2005, Manning et al. 2006, Barry 
and Giovannoni 2007). These previously characterized mutants 
were mostly selected from a production field or a breeding 
program as spontaneous variants (Giovannoni 2007). Tomato 
has six ethylene receptor genes in its genome (Klee 2004). 
Ethylene receptors control the fruit ripening process via the 
perception of ethylene. Nr is a well-characterized ethylene re- 
ceptor mutant that exhibits a strong ethylene-insensitive 
phenotype and inhibition of fruit ripening (Lanahan et al. 
1994, Yen et al. 1995). No new ethylene receptor mutant has 
been isolated in the past decade from tomato. 

In this study, the Micro-Tom TILLING platform was evalu- 
ated as a reverse genetics tool for the isolation of novel ethylene 



receptor mutant alleles. To attain this objective, six ethylene 
receptor genes {SIETR1-SIETR6) were screened by TILLING in 
our Micro-Tom mutant population. As a result, two Sletrl 
mutant alleles {Sletrl-1 and Sletr1-2) with reduced ethylene 
responses were identified. Delayed fruit ripening and prolonged 
fruit shelf life were observed in these mutants. We report that 
the Micro-Tom TILLING platform serves as a powerful tool for 
studying fruit biology and accelerates the isolation of important 
agronomic traits in tomato. The use of this tool as a new meth- 
odology is discussed. 



Results 

Mutation screening of a Micro-Tom EMS mutant 
collection by TILLING 

To apply TILLING methodology to the Micro-Tom mutant col- 
lection, DNA samples were prepared from 3,052 EMS- 
mutagenized M2 lines, generated in the presence of either 
0.5% EMS (2,180 lines) or 1.0% EMS (872 lines) (Watanabe 
et al. 2007, Saito et al. 2011) (Fig. 1). DNA superpools consist- 
ing of eight families per pool were subjected to a first screening. 
We used the CODDLE (Codons Optimized to Discover 
Deleterious Lesions) program for primer design to select opti- 
mal target amplicons, as performed previously (Dahmani- 
Mardas et al. 2010, Piron et al. 2010). All steps of the mutation 
screening were carried out by the procedure shown in Figs 1 
and 2. To estimate the mutation frequency of the Micro-Tom 
mutant collection and to assess its utility in the study of fruit 
biology, we screened for mutations in 10 genes involved in fruit 
ripening, softening and y-aminobutyric acid (GABA) metabol- 
ism. A total region of the tomato genome of approximately 
15.3kb was screened and 35 mutations were identified, 
among which 18 were from the 0.5% EMS population and 17 
were from the 1.0% population (Table 1). Tansition from G/C 
to A/T was observed in 91.4% of the mutations, known as the 
major nucleotide substitution induced by EMS (Greene et al. 
2003), while non-G/C to A/T substitution was also observed at 
low frequency (Supplementary Table S2). The mutation fre- 
quency was calculated as follows: the size in base pairs of the 
region screened was multiplied by the total number of lines 
screened per total number of identified mutations. The 0.5% 
EMS population had a mild mutation frequency of one muta- 
tion per 1,710 kb, whereas the frequency of the 1.0% EMS popu- 
lation was one mutation per 737 kb (Table 1). One to six 
mutated alleles were obtained per screened gene. Among the 
exonic mutations, the ratio of silent and missense mutations 
was 16.0 and 80.0%, respectively, and the ratio of nonsense 
mutations causing a stop codon was 4.0% (Table 2). 

Identification of novel Sletrl mutant alleles 
caused by EMS mutagenesis 

Mutations in the transmembrane region of ethylene receptor 
genes are known to affect the ethylene sensitivity in both 
Arabidopsis and tomato dramatically (Chang et al. 1993, 
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Fig. 1 Outline of Micro-Tonn TILLING platfornn. Seeds were nnutagenized with 0.5 or 1.0% EMS. To avoid ambiguities caused by chinnerisnn of 
mutant plants in the Mi generation, plants were self-fertilized. The M2 progeny of 10 seeds per line was used for screening. The DNA pool and M3 
seeds were collected from each M2 line (10 plants/line). An 8x DNA pool was used for the screening. PGR was performed using 5' end-labeled 
gene-specific primers to target the desired region, and heteroduplexes were formed by heating and cooling the PGR products. Recombinant 
SIEND01 nuclease was used to cleave base mismatches, and products that represent induced mutations were visualized by denaturing PAGE. 
Recombinant SIEND01 was extracted and purified from SIENDOI-overexpressing transgenic tomato plants. 



Wilkinson et ai. 1995). Therefore, by means of TILLING, we 
expect to screen effectively for and isolate various novel 
mutant alleles of ethylene receptor genes with weak and 
strong phenotypes and even for those with loss of function. 
We performed TILLING screening of six ethylene receptor genes 
in tomato, SIETR1-SIETR6, and found more than one allele for 
each gene (Table 1). A possible loss-of function mutant caused 
by a nonsense mutation (Q368stop) was found for SIETR5 
(Supplementary Table SI). Two ETR1 mutant alleles 
{Sletrl'1 and Sletrl-l) were found; each had an amino acid 
substitution in the predicted transmembrane region, c152t 
(P51L) for Sletrl'1 and t206a (V69D) for Sletrl-l (Fig. 3). The 
P51 L substitution in Sletrl-I corresponds to the same mutation 
position as Neuer ripe, where P36 was substituted to L in Nr 
(Wilkinson et al. 1995) (Supplementary Fig. S2), suggesting 
that the Sletrl-I mutation also affects receptor function and 
ethylene sensitivity. The hydrophobicity of the transmembrane 
region and the charge distribution in the flanking regions 
appear to be important determinants of the membrane orien- 
tation (Harley et al. 1998). The V69D substitution in Sletrl-l 
was located in the transmembrane-flanking region, as shown in 
the topological model of GmERSI proposed by Ma et al. (2006). 
The average hydrophobicity of Sletrl-l decreased from 0.176 in 
the wild-type protein to 0.166, as calculated by the SOSUI pro- 
gram (http://bp.nuap.nagoya-u.ac.jp/sosui/) (Hirokawa et al. 
1998), suggesting that this second allele of V69D would also 



affect the receptor function. Thus, we identified two 
Sletrl mutant alleles that possibly exhibited altered protein 
function. 

Sletrl mutant alleles resulted in ethylene 
insensitivity at different levels 

In this study, the ethylene-related phenotypes of the Sletrl 
mutant alleles {Sletrl-1 and Sletrl-l) were further character- 
ized. We tested whether the M3 generation of the mutant 
exhibited the triple response to exogenously applied ethyl- 
ene. Plants carrying the mutant alleles exhibited the 
ethylene-insensitive phenotype, whereas wild-type seedlings 
became swollen, and hypocotyl and root elongation were 
inhibited in response to ethylene treatment (Fig. 4A). 
Similar to Arabidopsis etr7-7, Sletrl-I seedlings were com- 
pletely ethylene insensitive (Ghang et al. 1993, Hua et al. 
1995), and failed to exhibit a reduction in hypocotyl and 
root length or an apical hook in response to ethylene treat- 
ment. The length of hypocotyls and roots were 58.7 ±1.9 
and 53.3 ±4.7 mm in the absence of ethylene and 57.7 ±0.5 
and 56.0 ±2.8 mm in the presence of ethylene (mean±SD), 
respectively. Sletrl-1 displayed a relatively weak ethylene 
insensitivity, and exhibited an intermediate level of hypo- 
cotyl length reduction compared with the wild type and 
Sletrl-I, from 50.7 ±2.5 to 31.7 ±2.4 mm in the Sletrl-l 
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Fig. 2 Work flow of TILLING in this study. TILLING was performed in this study as described in the figure. See Supplementary Fig. SI for an 
image of the results of TILLING screening. 



Table 1 Mutation frequency in the Micro-Tom EMS mutant population 

Gene Screened size (bp) No. of screened M2 lines No. of identified mutations'' Mutation frequency'' 







0.5 % EMS 


1.0 % EMS 


0.5 % EMS 


1.0 % EMS 


0.5 % EMS 


1.0 % EMS 


S/E7R7 


1,403 


2,180 


872 


4 


2 


1/765 kb 


1/612 kb 


S/E7R2 


1,185 


2,180 


872 


2 


0 


1/1,292 kb 




S/ETR3 (NR) 


945 


2,180 


872 


0 


2 




1/412 kb 


S/E7R4 


1,466 


2,180 


872 


0 


1 




1/1,278 kb 


S/E7R5 


1,447 


2,180 


872 


1 


4 


1/3,154 kb 


1/315 kb 


S/E7R6 


1,456 


2,180 


872 


0 


1 




1/1,269 kb 


SiSSADH 


1,572 


2,180 


872 


1 


2 


1/3,427 kb 


1/685 kb 


SIGABAT1 


3,537 


2,180 


872 


7 


4 


1/1,102 kb 


1/771 kb 


SIGABAT3 


1,007 


1,516 




2 




1/763 kb 




SIPL 


1,348 


664 


872 


1 


1 


1/865 kb 


1/1,175 kb 


Total/mean 


15,366 






18 


17 


1/1,710 kb 


1/737 kb 



THe number of identified mutation is shown based on TILLING gel screening. 
^ The mutation frequency is calculated as follows: (size in base pairs of the screened region) x (total number of lines screened)/(total number of identified mutations). 
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seedlings. The reduction of root length in Sletrl-l was simi- 
lar to that in Sletrl-I in the presence of lOp.p.m. ethylene, 
showing only a slight reduction, from 573 ±2.1 to 
46.3 ±2.4 mm in the Sletrl-l roots (Fig. 4A, B). The epinas- 
tic response exhibited a similar trend to the triple response. 
In the presence of SOp.p.m. ethylene, the petiole exhibited 
downward curvature in the wild type, but showed a weaker 
or no response in Sletrl-I and Sletrl-l (Fig. 4C). 

Sletrl mutant alleles exhibited dominant 
inheritance 

Constitutive overexpression of mutated ethylene receptor 
genes and genetic analyses in Arabidopsis and tomato have 
revealed that ethylene receptor mutants inherit the ethylene- 
insensitive phenotype in a dominant-negative manner (Chang 
et al. 1993, Lanahan et al. 1994, Wilkinson et al. 1995). To in- 
vestigate the inheritance pattern of the obtained Sletrl mu- 
tants, we crossed the two mutants with wild-type Micro-Tom 
and Ailsa Craig. Sletrl-1 and Sletr1-2 F2 seedlings were treated 

Table 2 The type of amino acid substitution in the exonic region 



Effect No. of Percentage (%) 

mutations 



Silent mutation 


4 


16.0 


Missense mutation 


20 


80.0 


Nonsense mutation 


1 


4.0 


Total 


25 


100 



with 10 or 5 p.p.m. ethylene. The triple response was assayed for 
each BCtSi population (S/etr7-7 x Micro-Tom BC1S1, Sletrl- 
2xMicro-Tom BC1S1) and the result was scored as ethylene 
sensitive (5), intermediate (9) or insensitive (6) for the 
Sletrl-1 BC1S1 population (x^ = 0.3, P > 0.8), and ethylene sen- 
sitive (10), intermediate (18) or insensitive (10) for the Sletrl-2 
BC1S1 population (x^ = 0.237, P>0.8). The F2 population 
(S/etr7-7x Ailsa Craig F2, S/etr7-2x Ailsa Craig F2) was scored 
as ethylene sensitive (10) or insensitive (33) for the Sletrl-1 F2 
population (x^ = 0.07, P > 0.75), and ethylene sensitive (26) or 
insensitive (63) for the Sletrl-2 F2 population (x^ = 0.843, 
P > 0.3) (Table 3). The F2 populations were scored only as 
ethylene sensitive and insensitive, because they exhibited an 
indistinguishable intermediate phenotype of different seedling 
lengths due to the effect of the Micro-Tom or Ailsa Craig gen- 
etic background. The segregation ratio of both of the BC1S1 and 
F2 populations in Sletrl-1 and Sletrl-2 was nearly 1:2:1 in the 
BC1S1 and 3:1 in the F2, respectively, suggesting that the pro- 
geny inherited both alleles as a single dominant trait. The in- 
heritance pattern of the Sletrl mutant alleles was consistent 
with a previous report of the ethylene receptor mutant Nr, 
which displayed a semi-dominant triple response phenotype 
(Lanahan et al. 1994). 

Sletrl mutant alleles resulted in delayed petal 
abscission and fruit ripening 

Both Sletrl mutants showed delayed petal abscission (Fig. 5B). 
This phenotype is similar to that of the Nr mutant and RNAi 



hisietrl-l) 

MGSLLRMNRLLSSIVESCNCIIDPQLPADDLLMKYQY ISDFFIALAYFSIPVELIYFVKK 
D{Sletrl-2) 

SAVFPYRWVLVQFGAFIVLCGATHLINLWTFNMHTRNVAIVMTTAKA LTALVSCITALML 

W(2442) 

VHIIPDLL SVKTRELFLKKKAAQLDREMGIIRTQEETGRHVRMLTHEIRSTLDRHTILKT 

1(2162) 

TLVELGRTLALEECALWMPTRTGLELQLSYTLRHQNPVGLTVPIQLPVINQVFGTNHWK 

l"(4i"60) Fiilio) 
ISPNSPVARLRPAGKYMPGEVVAVRVPLLHLSNFQINDWPELSTKRYALMVLMLPSDSAR 



QWHVHELELVEVVADQVAVALSHAAILEESMRARDLLMEQNVALDLARREAEMAVRARND 



FLAVMNHEMRTPMHAIIALSSLLQETDLTPEQRLMVETILKSSNLLATLINDVLDLSRLE 
DGSLQLDIGTFNLHALFREVHSLIKPIASVKKLFVTLSLSSDLPEYVIGDEKRLMQILLN 
WGNAVKFSKEGNVSISAFVAKSDSLRDPRAPEFFAVPSENHFYLRVQIKDTGIGITPQD 
I PNLFSKFTQSQALATTNSGGTGLGLAI CKRFVNLMEGH I WI ESEGLGKGSTAI F I I KLG 
IPGRANESKLPFVTKLPANHTQMSFQGLKVLVMDENGVSRMVTKGLLTHLGCDVTTVGSR 
DECLRVVTHEHKWIMDVSMQGIDCYEVAWIHERFGKRHGRPLIVALTGNTDRVTKENC 

MRVGMDGVILKPVSVYKMRSVLSELLEHGVVLES 

Fig. 3 The amino acid sequence of the tomato ethylene receptor (SIETRI) and sequence variants identified by TILLING. Vertical lines and a 
dotted line indicate the transmembrane region and GAF domain, respectively. 
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Fig. 4 The ethylene response of Sletrl mutant alleles. (A) Triple response phenotype of the wild type (WT), Sletrhl and Sletrl-I seedlings. Seeds 
were surface sterilized and sown on 1/2 MS medium in the absence (-) or presence (+) of 10 p.p.m. C2H4 and incubated at 25°C in the dark for 
7d. Bar = 1cm. (B) Quantification of ethylene-induced inhibition of root and hypocotyl growth. Vertical bars represent the SD (n = 3). 
(C) Epinastic response phenotype of 4-week-old wild-type (WT), Sletrl-l and Sletrl-I plants. Plants were treated with 50 p.p.m. ethylene at 
25°C for 20 h. Bar =5 cm. 



Table 3 Inheritance pattern of Sletrl mutant alleles 



Population" F2 and BC1S1 Inheritance 

segregation*' value^ pattern^ 
(mutant : WT) 



S/etr7-7xAilsa Craig F2 


33 


10 


0.07 


Single 


dominant 


S/etr7-7x Micro-Tom BCiSt 


15 


5 


0 


Single 


dominant 


S/etr7-2xAilsa Craig F2 


63 


26 


0.843 


Single 


dominant 


S/etr7-2x Micro-Tom BCqSi 


28 


10 


0.035 


Single 


dominant 



Sletrl mutant alleles were crossed with wild-type (WT) Micro-Tom and Ailsa 
Craig. 

The number of progeny exhibiting WT (ethylene-sensitive) or mutant 
(ethylene-insensitive) phenotypes in the F2 and BCtSt populations is shown. 
The number of progeny exhibiting intermediate ethylene insensitivity in BCiSq 
is shown in the text. 
The values were calculated for the F2 and BCtSt populations. 

Inheritance pattern was estimated based on the value. The values were 
significant at the level of 5%. 

(RNA interference) transgenic lines of tomato ETHYLENE 
INSENSITIVE LIKE (EIL) (Lanahan et al. 1994, Yokotani et al. 
2009), which are associated with inhibition of the ethylene sig- 
naling pathway. Flowers with fully opened petals were defined 



as entering anthesis and tagged as 0 days after flowering (DAF). 
Wild-type Micro-Tom flowers began to wilt at 3 DAF, whereas 
flowers of S/etr7 mutants did not wilt even at 15 DAF, except for 
flowers of heterozygous Sletrl-l which wilted at 5 DAF 
(Fig. 5B). Compared with Sletrl-l, Sletrl-I exhibited a more 
severe delayed abscission phenotype, and petals were still at- 
tached to the fruits even at 60 DAF (Fig. 5A). Delayed fruit 
ripening was observed in Sletrl-I, and the fruit turned fully yel- 
low or orange 7-10 d later than did the wild type (Fig. 5A). 
At 60 DAF, most Sletrl-1 fruits did not turn red, while some 
changed to an intermediate color between orange and red 
(Supplementary Fig. S4). No difference was observed between 
Sletrl-l and the wild type regarding the ripening phenotype 
(Fig. 5A). These results indicate that the degree of petal abscis- 
sion and fruit ripening were associated with the level of ethyl- 
ene sensitivity. 

To confirm the linkage between the mutant alleles and 
phenotypes, we evaluated the abscission and ripening pheno- 
type of each BC1S1 population. Genotyping of Sletrl mutant 
alleles was carried out by cleaved amplified polymorphic se- 
quence (CAPS) analysis (see the Materials and Methods). 
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Petal abscission and the ripening phenotype in homozygous 
and heterozygous plants were compared with those in the 
wild type. Heterozygous plants exhibited an intermediate 
phenotype between wild-type and homozygous plants. Petals 
of heterozygous Sletrl-l plants began to wilt at 5 DAF, and 
almost separated from fruits at 7 DAF and 10 DAF (Fig. 5B). 
At the late ripening stage (50-60 DAF), heterozygous Sletrl-I 
fruits became red (Fig. 5C), indicating that the petal abscission 
and ripening phenotypes were semi-dominant in mature 
plants. The Sletrl-I and Sletrl-l mutations co-segregated with 
delayed petal abscission and ripening in each F2 population 
(data not shown). 

Fruit shelf life was prolonged in Sletrl-l 

Fruit shelf life is an important property in determining the 
market value of tomato. We hypothesized that weak ethylene 
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Fig. 5 Fruit ripening and petal abscission phenotype of Sletrl mutant 
alleles. Comparison of fruit ripening in the wild type (WT) and Sletrl 
mutant alleles. Fruits were harvested at different maturities. (B) Petal 
abscission phenotype in homozygous and heterozygous Sletrl mu- 
tants. (C) Fruit ripening phenotype in the homozygous and hetero- 
zygous Sletrl'1 mutant. Bar = 2 cm. Homozygous and heterozygous 
plants were selected by CAPS analysis (shown in the Materials and 
Methods). 



insensitivity might increase fruit shelf life, without visibly affect- 
ing the ripening color. The Sletrl-l mutant was an ideal plant to 
test this hypothesis. Wild-type and Sletrl-l fruits of the same 
stage (45 DAF) were harvested and stored in a sealed chamber 
at 25°C for 60 d. In wild-type fruits, black spots began to arise on 
the fruit surface at 20-25 d after harvest, whereas the Sletrl-l 
fruits did not show any spots (data not shown). The wild-type 
fruits deteriorated 60 d after harvest, whereas in the Sletrl-l 
fruits the surface remained intact (Fig. 6). This observation 
shows that fruit shelf life is dramatically prolonged in Sletrl-l 
and that deterioration of post-harvest fruit is delayed. 



Discussion 

The dwarf cultivar Micro-Tom has recently been identified as a 
model tomato (Matsukura et al. 2008). Genome information of 
Micro-Tom, such as full-length cDNA and bacterial artificial 
chromosome (BAC) sequences, has been accumulated by the 
Kazusa DNA Research Institute as part of the National 
BioResource Project (NBRP, MEXT) (Aoki et al. 2010). We pre- 
viously developed a Micro-Tom mutant population and a data- 
base named 'TOMATOMA', which provides information about 
visible mutant phenotypes. 

In this study, we generated a Micro-Tom TILLING platform 
from 3,052 EMS-mutagenized lines. We also deployed the 
LI-COR/EN D01 system for detection of mutations, as a 
low-cost and high-sensitivity method (Fig. 1) (Triques et al. 
2007, Triques et al. 2008, Piron et al. 2010). To be cost-effective, 
we extracted endonuclease I for digesting mismatch DNA from 
SIENDOI-overexpressing transgenic tomato plants (see the 
Materials and Methods). The mismatch DNA recognition and 
cleavage activity of SIEND01 was evaluated by using different 
combinations of artificially designed mismatches. We could see 
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Fig. 6 Fruit shelf life of the Sletrl-l mutant. Fruits were harvested 45 d 
after flowering and stored at 25°C for 60 d in a sealed chamber. The 
appearance of Sletrl-l fruits was compared with that of the wild type 
at 60 days after harvest. Bar = 2 cm. 
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no bias in DNA mismatch types recognized and cleaved by 
SIEND01 (Supplementary Fig. S5A), indicating a similar ten- 
dency to the Arabidopsis END01 (Triques et al. 2007). The 
transgenic tomatoes expressing SIEND01 stably provided the 
enzyme with a simple extraction procedure; we used this 
system in the TILLING analyses. Screening of 10 genes by tar- 
geting an approximately 15.3 kb region of the tomato genome 
has enabled us to estimate the mutation frequency to be one 
mutation per 1,710 kb for the 0.5% EMS population, and one 
mutation per 737 kb for the 1.0% EMS population. A correlation 
was observed between the concentration of EMS applied and 
the mutation frequency obtained, where the frequency was 
2.32 times higher in the 1.0% population than in the 0.5% 
population. This result is consistent with the previously re- 
ported Red Setter mutant population (Minoia et al. 2010). 
The mutation frequency obtained for our 1.0% EMS population 
was lower than that of the Red Setter population (one per 
322 kb; Minoia et al. 2010), but was similar to that of the 
TPAADASU population (one per 737 kb; Gady et al. 2009) 
and the M82 population (one per 574 kb; Piron et al. 2010). 
Thus, we concluded that the 1.0% populations are more effi- 
cient for screening mutants by applying TILLING using our 
platform. We intend to increase the size of the 1.0% population, 
and thereby increase the chance of isolating mutants of target 
genes. 

Given that one would expect a higher concentration of 
EMS to result in a higher mutation frequency and to facili- 
tate the identification of desired mutant alleles even in a 
small-scale population, more backcrossing might be neces- 
sary to characterize the specific mutant phenotypes due to 
the presence of mutations that might affect unpredictable 
plant growth. In fact, the Micro-Tom mutants recovered 
from the populations with higher EMS concentrations 
tend to show stronger infertility, most probably due to 
the multiple mutations that could affect reproductive de- 
velopment (Saito et al. 2011). Thus we suggest that the 
mutants isolated from our TILLING platform should be sub- 
ject to backcrossing several times to remove undesirable 
mutations before mutant analysis. 

Recently, a practical approach using identified mutants 
has been made by taking advantage of the TILLING plat- 
forms developed in many crop species. For example, inter- 
esting mutants were identified in tomato, namely the 
SlelF4E1 G1485A mutant, which showed potyvirus resistance, 
and in melon, namely the CmACO! G194D mutant, which 
produced fruit with an enhanced shelf life (Dahmani-Mardas 
et al. 2010, Piron et al. 2010). These reports indicated that 
TILLING is a useful and could be a practical tool for crop 
improvement. Our study also provided evidence that the 
Micro-Tom TILLING platform can be used as a valuable 
tool for studying fruit biology and isolating genetic material 
that harbors important agronomic traits in tomato. 
Currently, >2,000 new lines have been added to the 
Micro-Tom TILLING population to improve the quality of 
the platform. This population will become available in the 



near future and will improve the possibility of identifying a 
greater number of interesting mutants. 

As a result of SIETR1 screening, we identified two mutant 
alleles {Sletrl-I and Sletrl-l) that exhibited reduced ethylene 
sensitivity (Fig. 4). These mutations were located in different 
domains of the transmembrane region, i.e. P51L (Sletrl-I) was 
in the first and V69D (Sletrl-l) in the second domain (Fig. 3). 
The P51L substitution of Sletrl-1 corresponds to the same 
amino acid substitution of P36L in Nr and Arabidopsis etr2-1 
(Supplementary Fig. S2; Sakai et al. 1998). This proline residue 
in the first domain is highly conserved amongst known and/or 
putative ethylene-binding proteins (Wang et al. 2006). Nr is 
reported to show an ethylene-insensitive phenotype in the 
seedling triple response assay. In the dose response assay of 
exogenously applied ethylene, Nr responded to a low concen- 
tration of ethylene (1 p.p.m.) (Lanahan et al. 1994). In contrast, 
Sletrl-1 showed no response to 10 p.p.m. ethylene, indicating 
that this mutant is less sensitive to ethylene than is Nr. 
The ethylene insensitivity of Sletrl-1 resembled that of 
Arabidopsis etr1-1 (Atetr1-1), which is completely insensitive 
to ethylene. The relationship between amino acid substitu- 
tion and ethylene binding capacity was investigated in 
Arabidopsis, and it was found that the binding capacity of 
P36L-mutated AtETRI was <5% that of the wild-type protein 
(Wang et al. 2006). This result supports the hypothesis that 
P51L-mutated SIETR1 is incapable of binding to ethylene. 
The second affected allele, V69, was located near the ethy- 
lene-binding domain, where the effect on ethylene binding 
was previously reported (Rodriguez et al. 1999, Wang et al. 
2006). The hydrophobicity of the transmembrane region ap- 
pears to be an important determinant of the membrane orien- 
tation (Harley et al. 1998). The calculated hydrophobicity 
was reduced as a result of the V69D substitution in the 
Sletr1-2 protein. However, no experimental evidence has been 
reported to date that demonstrates the importance of this 
residue in ethylene binding activity. The fact that Sletr1-2 ex- 
hibited an ethylene-insensitive phenotype suggests that the 
V69D substitution of Sletr1-2 may also affect ethylene binding 
activity. 

Homozygous Sletrl-1 and Nr showed similar fruit ripening 
phenotypes, both being impaired in ripening at 60 DAP 
(Fig 5A). The ripening phenotypes of the heterozygous lines 
differed from those of the homozygous lines. Whereas the fruit 
of heterozygous Sletrl-1 finally turned red (Fig. 5C), the ripen- 
ing of Nr fruit remained incomplete with respect to color 
(Lanahan et al. 1994). Since we were concerned about the 
effect of the different genetic backgrounds of Sletrl (Micro- 
Tom) and Nr (Ailsa Craig), we evaluated the petal abscission 
phenotype in Nr (Nr/Nr) (cv. Rutgers) x wild type (nr/nr) 
(cv. Micro-Tom) BCi plants (Nr-MT BCi). We used Nr in the 
Rutgers background for this experiment, since Nr in the Ailsa 
Craig background was unavailable at the time. Heterozygous 
Sletrl-1 exhibited impaired petal abscission, whereas the effect 
was weaker for Nr-MT BCi (data not shown), indicating that 
the ripening and petal abscission phenotypes of heterozygous 
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Sletrl'1 were opposite to those of Nr-hAJ BCi. To test if we 
could obtain a similar result with Sletrl-I in a different genetic 
background we generated Ft heterozygous lines of Sletrl-I 
crossed with Ailsa Craig. These lines bore red fruits and ex- 
hibited an apparent delay in petal abscission (data not 
shown). Whereas Sletrl-I shows a stronger degree of ethylene 
insensitivity than Nr in the seedling, the opposite is true during 
the fruit ripening process. NR gene expression is up-regulated at 
the breaker stage together with the initiation of autocatalytic 
ethylene biosynthesis at the onset of ripening, whereas SIETR1 
gene expression stayed at a constant level during fruit develop- 
ment (Lashbrook et al. 1998). The gene expression pattern of 
SIETR1 and NR during fruit development (Lashbrook et al. 
1998) also supports this idea that NR is important in the fruit 
ripening process. Together with previous reports, our data sug- 
gest that SIETR1 and Nr are functionally distinct in petal abscis- 
sion and fruit ripening. Although no clear ripening phenotype 
was observed for Sletr1-2 (Fig. 5A), deterioration of post- 
harvest fruit was prolonged compared with that in the wild 
type (Fig. 6), indicating that the weak ethylene insensitivity 
affected an aspect of fruit ripening. Improvement of post- 
harvest fruit shelf life is an important agronomic trait in tomato 
breeding. Modifying the key regulatory factors of fruit ripening 
and softening may prolong shelf life and provides massive com- 
mercial and economical benefit. Ethylene controls the fruit 
ripening process of climacteric fruits via the ethylene biosyn- 
thetic and signaling pathway (Barry and Giovannoni 2007). 
While several ripening mutants, such as rin, nor and Nr, have 
been identified (Barry and Giovannoni 2007), only rin has been 
used in practical breeding applications (Kitagawa et al. 2005). 
Although Nr produces fruit with a good shelf life, it has not 
been used in breeding programs due to its incomplete ripening 
phenotype even in the heterozygous state (Nr/nr), which re- 
sults in insufficient red coloration (Lanahan et al. 1994). 
Furthermore, there are traits associated with Nr that are un- 
desirable, such as increased susceptibility to some pathogens 
(Francia et al. 2007, Kavroulakis et al. 2007, Cantu et al. 2009). 
Mutants with drastically altered phenotypes may be useful ma- 
terial for basic studies, but are often not suitable for use in 
breeding. In this regard, since Nr has a more pronounced im- 
pairment of ripening than does Sletrl, the Sletrl alleles are more 
suitable for use in breeding programs. The identification of two 
novel Sletrl mutant alleles that are distinguished by the level of 
ethylene sensitivity and the characterization of their associated 
phenotypes could provide insight into the ethylene-mediated 
fruit ripening mechanism in tomato. 

Materials and Methods_^^^^^^^^^^^^| 

Plant material 

Tomato (S. lycopersicum) cv. Micro-Tom was used. The 
TILLING population was developed as previously described 
(Watanabe et al. 2007, Saito et al. 2011). EMS-mutagenized 
plants were grown in a greenhouse. 



DNA extraction and pooling 

For each M2 line, genomic DNA was isolated from 10 individual 
plants. Leaves were collected in 2 ml micro test tubes, frozen in 
liquid nitrogen, and homogenized with a pestle. Genomic DNA 
was extracted using a Maxwell16 Tissue DNA Purification Kit 
(Promega). To generate DNA superpools in 96-well format, 
DNA from eight lines was mixed in one well. The concentration 
of DNA was normalized by separating DNA samples on an 
agarose gel by electrophoresis and comparing the intensity of 
the ethidium bromide-stained bands with that of a ADNA ref- 
erence (New England Biolabs). 

PCR amplification 

PGR primers were designed using the PRIMER3 program (Rozen 
et al. 2000) after selecting an optimal region to screen using the 
CODDLE program (Till et al. 2003) except for SISSADH and 
SIGABAT1 for which entire coding regions were targeted. PCR 
amplification was carried out by two methods, either nested 
PCR was performed using universal primers or gene-specific 
primers were used. In the gene-specific primer system, PCR 
was carried out by direct PCR using IRD700 and IRD800-labeled 
primers; in the universal primer system, non-labeled gene- 
specific primers attached to the T7 (CGCGTAATACGACTCA 
CTATAG) or SP6 (CATACGATTTAGGTGACACTATAG) se- 
quence at the 5^ end were used in the first round of PCR. 
PGRs were performed in a total volume of 15 |il, with 5-10 ng 
of superpool DNA. DNA was mixed with 1.5 |il of 5x Green 
Gotaq Flexi buffer, 2.0 |il of 25 mM MgCl2, 1.2 |il of 2.5 mM 
dNTP, 0.3 jil of 10 |iM primers and 0.075 |il of 5 U jil"^ Gotaq 
Hotstart DNA polymerase or Gotaq Flexi DNA Polymerase 
(Promega). Then, 0.5 |il from the first PCR was used as a 
template in the second nested PCR using IRD700- and 
IRD800-labeled T7 and SP6 primers (biomers.net). After 
confirming successful PCR amplification by electrophoresis, 
3-7 |il of PCR products was used in the subsequent TILLING 
screening. 

Preparation of tomato END01 nuclease and 
detection of mismatch DNA 

Coding sequences of SIEND01 (accession No. AB667996) were 
cloned into the pGWB8 binary vector (Nakagawa et al. 2007). 
Recombinant SIEND01 tagged with six histidine residues 
was stably expressed in tomato cv. Micro-Tom. SIEND01 pro- 
tein was extracted from transgenic tomato leaves. We basically 
followed a previously described protein purification method 
(Triques et al. 2007), with some modifications. A 2g aliquot 
of tomato leaves was homogenized with a mortar and pestle 
in 10 ml of buffer [0.1 M Tris-HCI (pH 8.0), 1 mM phenyl- 
methylsulfonyl fluoride (PMSF), 0.125 mM 2-mercaptoethanol 
and 10% glycerol]. The extract was cleared by centrifugation 
at 16,000xg for 20min. The supernatant was used for fur- 
ther purification. Recombinant SIEND01 was purified using a 
Histrap HP Column (GE Healthcare). The eluted fraction 
was dialyzed as previously described (Triques et al. 2007). The 
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enzyme activity for each batch of SIEND01 was determined 
using a dilution series of control samples (Supplementary 
Fig. S5B). PCR products were mixed with sterilized water 
up to a total volume of 10 |il and subjected to SIEND01 
digestion as previously described (Triques et al. 2007). 
TILLING screening was performed using an LI-COR DNA ana- 
lyzer as previously described (Till et al. 2006) and as shown 
in Fig. 1. 

Selection of TILLING mutant alleles in bulked M3 
line populations 

To select mutant individuals, a CAPS marker was designed for 
each identified allele to distinguish the mutant allele from that 
of the wild type. Appropriate restriction enzymes were selected 
from those listed by the PARSESNP tool (Taylor and Greene. 
2003), as shown in Supplementary Fig. S3. If no restriction site 
was found at the position of mutation, the target region was 
amplified with TILLING primers used in the screening, and the 
obtained PCR products (400-500 ng) were digested with 
SIEND01. Digested fragments were visualized by electrophor- 
esing in a 1.5-2.0% standard agarose gel followed by SYBR Safe 
DNA gel staining (Invitrogen). 

Seedling triple response assay and epinastic 
response assay 

Seeds were sterilized with 10% commercial bleach including a 
detergent (Kitchen Haiter, Kao) for 20 min and then rinsed with 
sterilized water three times for 5 min each. The seeds were 
germinated in a 50 ml glass bottle containing 10 ml of 1/2 MS 
medium (Murashige and Skoog 1962). Ethylene was added to 
the bottles sealed with silicon rubber at the designed concen- 
trations (5 or 10 p.p.m.) and seedlings were grown for 7 d in the 
dark at 25°C. The epinasitic response assay was performed 
using 4'Week'Old plants, ethylene was injected into the sealed 
chambers at the designed concentration (50 p.p.m.), and the 
plants were incubated in growth chambers for 20 h in the light 
at 25°C. 

Characterization of mutant phenotypes 

Seeds of wild-type Micro-Tom and BC^S^ and BCtS2 homozy- 
gous seeds of Sletrl-I and Sletrl-l were sown on wet filter paper 
and placed under continuous light for 2-3 d at 25°C to stimu- 
late seed germination. Germinated seeds were then trans- 
planted into soil and grown in a photoperiod of 16 h light at 
25°C. The date of flowering (0 DAF) was tagged for the time 
course observation of flower abscission, fruit ripening and 
evaluation of fruit shelf life at the times indicated in Figs. 5 
and 6. All plants were grown under the same conditions as 
described above. 

Genotyping of Sletrl mutant alleles 

Homozygous and heterozygous Sletrl plants were selected from 
BCtSt populations as follows. CAPS analysis was used to distin- 
guish the Sletrl'1 or Sletrl-l allele from the wild-type SIETR1 



allele. PCR amplification for the detection of each mutation was 
performed with a SIETR1-CAPS forward primer (5^-gtataaaagg 
agttggggcaaag-30 and SIETR1-CAPS reverse primer (5^-atcagg 
aatgatgtggacaagc-30 to yield a 730 bp PCR product. PCR prod- 
ucts were digested with Xspl (TAKARA) and A/lbol (New 
England Biolabs) to detect Sletrl-I and Sletrl-l alleles, respect- 
ively. The Xspl-digested Sletrl-1 product consisted of a 503 and 
227 bp fragment. Xspl digestion of SIETR1 /Sletrl-1 heterozygous 
plants resulted in three fragments, the 730 bp fragment 
(non-digested wild-type allele) and a 503 and 227 bp fragment 
(digested Sletrl-1 allele). /VI bo I -digested Sletr1-2 consisted of a 
159, 397 and 174 bp fragment. Mbo\ digestion of SIETR1 /Sletrl -2 
heterozygous plants resulted in four fragments, of 159 and 
571 bp (partially digested wild-type allele) and of 397 and 
174 bp (digested Sletrl-2 allele). 
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